Low-background gamma counting at the Kimballton Underground Research Facility 
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The next generation of low-background physics experiments will require the use of materials with 
unprecedented radio-purity. A gamma-counting facility at the Kimballton Underground Research 
Facility (KURF) has been commissioned to perform initial screening of materials for radioactivity 
primarily from nuclides in the 238 U and 232 Th decay chains, 40 K and cosmic-ray induced isotopes. 
The facility consists of two commercial low-background high purity germanium (HPGe) detectors. 
A continuum background reduction better than a factor of 10 was achieved by going underground. 
This paper describes the facility, detector systems, analysis techniques and selected assay results. 

PACS numbers: 07.85.Fv, 29.40.-n, 07.85. Nc, 29.30.Kv 
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I. INTRODUCTION 

We have commissioned a low-background gamma- 
counting facility at the Kimballton Underground Re- 
search Facility (KURF). KURF is located at Lhoist 
North America's Kimballton mine in Ripplemead, Vir- 
ginia. The experimental hall is located on the mine's 14 th 
level at a depth of 1450 m.w.e (meters of water equivalent 
shielding) . The overburden consists of 520 m of dolomite, 
limestone and other sedimentary rock. Experiments are 
housed in a 30 m x 11 m laboratory building that was 
completed in October 2007 (Fig. [TJ. The laboratory's 
general infrastructure is maintained primarily by collab- 
orators at Virginia Polytechnic Institute and State Uni- 
versity. KURF has an office, air filtration, power, water, 
phone and ethernet. KURF also has the advantage of 
drive in access, making it simple to transport personnel 
and equipment to the experimental hall. Liquid nitro- 
gen (LN 2 ) used in the laboratory is stored in a 2.4 m 3 
portable dewar that can be transported to the surface to 
be refilled as necessary. Radon levels in the laboratory 
have been found to vary from 37 Bq/m 3 in the winter to 
122 Bq/m 3 in the summer. The detectors used in this 
paper are housed in sealed modified shipping containers 
(MSCs) within the laboratory building (Fig. [I). 



Series coaxial detector [l|; the high voltage filter and 
preamplifier are removed from the line of sight of the 
crystal to reduce backgrounds from radioactive contami- 
nants that may exist in detector components. The sample 
cavity is cylindrical, 41 cm (height) x 28 cm (diameter). 
Further specifications for VT-1 are in Table El A VT-1 
background spectrum taken on the surface and a back- 
ground spectrum taken underground is shown in Fig. O 

The second detector, named "Melissa," is a Canberra 
LB (low-background) coaxial detector 0). Melissa is in 
a vertical orientation with a dipstick style cryostat. The 
preamplifier is also removed from the cryostat, allowing 
for shielding to be placed in between the preamplifier and 
cryostat. The sample cavity is 38 cm x 38 cm x 38 cm. 
Further specifications for Melissa are in Table [J 

A Melissa background spectrum is shown in Fig. [3] 
No active shielding is currently used for cither detector. 
The sample cavities of both detectors arc continuously 
purged with LN 2 boil-off to flush out radon. A reduction 
of ^80% in the activity of radon daughters was observed 
after introducing the LN 2 boil-off, however activity due 
to radon remains as a dominant source of background. 
The shield design of the detectors has made it difficult to 
design a hermetically sealed radon exclusion system. 



II. DETECTORS &i SHIELDING 

The counting facility consists of two high purity ger- 
manium (HPGe) detectors specifically designed for low- 
background assay work. The first detector, named "VT- 
1" , is a commercial ORTEC LLB (very low-background) 
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Integral count rates for these background spectra 
(Figs. [H and [3J arc shown in Table ITT1 While direct com- 
parisons with other underground gamma-counting facili- 
ties are difficult, integral count rates for other selected fa- 
cilities can be found in [J-Q ■ Although the integral count 
rate in the 40-2700 keV region is higher than for other 
detectors at a similar depth, most of our background is at 
low energy. An overview of low-radioactivity background 
techniques and a comparison of low-level counting meth- 
ods can be found in [7], [8[ and references therein. 




FIG. 1: KURF before and after construction of the laboratory building. Modified shipping containers (MSCs) are 

shown in the right panel (Color online). 



TABLE I: Detector specifications. 





Melissa 


VT-I 


Manufacturer 


Canberra 


ORTEC 


Relative Efficiency 


50% 


35% 


Performance 






FWHM at 1.33 MeV (keV) 


1.70 


1.80 


Threshold (keV) 


20 


20 



Shield 



Crystal (coaxial) 



Cryostat 



Lead thickness (cm) 
Oxygen-free high conductivity 
(OFHC) copper thickness (cm) 

Mass (kg) 

Length (mm) 

Diameter (mm) 

Hole diameter (mm) 

Hole depth (mm) 

Outer electrode thickness (mm) 

Inner electrode thickness (/xm) 

End cap diameter (mm) 
End cap thickness (mm) 
End cap to crystal (mm) 
End cap material 
IR window material 



15.2 



2.54 



10. 1 



0.3 



1.1 


0.956 


64.5 


75.7 


65 


55.8 


7.5 


9.1 


50 


63.2 


1.06 


0.7 


0.3 


0.3 


82 


70 


4.2 


1.3 


5 


4 


igh purity Al 


Mg 



Al Mylar iM , Kapton Al Mylar 



TM 



TABLE II: Integral count rates for Melissa and VT-1 
in [10 3 counts/day]. 

Melissa VT-1 VT-1 (surface) 



40-2700 keV 
40-1000 keV 
1000-2700 keV 



7.8 
7.5 
0.27 



7.6 
6.5 
0.71 



68 
16 



III. SAMPLE PREPARATION 

To achieve the required assay sensitivity for next gen- 
eration low-background experiments, unwanted radionu- 
clides must be removed from sample surfaces to prepare 
for gamma-counting. Samples are prepared for assay in 
a clcanroom environment using ultra-pure reagents and 
clean plastics that have been screened for radioactivity. 
Depending on the sample material and required detection 
limit, a variety of methods can be used to treat sample 
surfaces. These methods include acid leaching for plas- 
tics, acid etching for metals or cleaning with ultra-pure 
solvents. After cleaning, samples arc bagged in nylon to 




Energy [keV] 

FIG. 2: VT-1 surface (black/top) and underground 
(red/bottom) comparison. Both spectra were taken 

with the lead shield in place (Color online). The 

underground spectrum was taken with a LN2 boil-off 

purge in place 
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FIG. 3: Background spectrum for the Melissa 

detector. The spectrum was taken with lead shielding 

and a LN2 boil-off purge in place (Color online) . 



dards [T(|. This process is complicated and time con- 
suming and is of limited accuracy for complex geometries. 
Another method is to use a point source calibration at 
representative points determined by Monte Carlo simu- 
lation. The efficiency curve generated is then corrected 
for absorption by a sample matrix and sample container 

B 

At KURF, a detailed Monte Carlo simulation for each 
sample is done using MaGe [l2|. MaGe is a Geant4- 
based [13|, [l4| simulation package maintained and devel- 
oped by a joint group of the Majorana [15l - tr7j and 
Gerda [18| collaborations. Once a detailed sample ge- 
ometry has been coded into the simulation, it is uniformly 
doped with isotopes of interest from the 238 U, 232 Th, and 
40 K decay chains, along with any other isotopes that may 
be present in the sample, e.g. 60 Co. The primary 7-rays 
from these decays are tracked from the emission at the 
source to absorption in the detector active region. By 
using a pure Monte Carlo simulation to determine the 
FEP detection efficiency, self-attenuation in the sample 
is accounted for and there are no limitations on source 
or detector configurations. 

The simulated spectrum is then convolved with the fi- 
nite energy resolution of the detector. The energy resolu- 
tion for both detectors has been measured as a function of 
7-ray energy from 303-1836 keV using radioactive point 
sources. 

The peak area in the simulated spectrum is deter- 
mined by fitting the peak of interest with a Gaussian 
and subtracting a linear background. The FEP detec- 
tion efficiency, including branching ratios, can then be 
determined from 



peak area 



number of events simulated 



B. Monte Carlo validation 



(1) 



prevent any recontamination of surfaces. 



IV. ANALYSIS TECHNIQUES 

A. Monte Carlo simulations and efficiency 
calculations 

The full energy peak (FEP) detection efficiency, de- 
fined as the ratio of the number of events detected in the 
gamma peak to the number of events emitted from the 
source for a particular radioactive isotope in a sample, 
depends on many factors, including the crystal, cryo- 
stat, shielding and source geometries. There are cur- 
rently several methods to determine the FEP detection 
efficiency. One is to use analytical calculations Q, how- 
ever this technique is limited to simple geometries and 
requires complex calculations. In some cases, a physical 
model of the sample can be created using known stan- 



It is our goal to know the efficiency of any volume 
sample to < 10% and attribute < 10% systematic un- 
certainty to these efficiencies. In order to validate the 
Monte Carlo FEP detection efficiency calculations, point 
sources of well-known activity were used to baseline the 
simulation using a method similar to [19j . Once the de- 
tector geometries are well understood, the Monte Carlo 
can be used to simulate all relevant physical processes 
and accurately determine peak efficiencies. 

Experimental data were taken for each detector using 
point sources emitting 7-rays with energies ranging from 
303-1836 keV. The activity of each source was known to 
within 3%. To understand the effects of source place- 
ment relative to the crystal, the sources were placed in 
the locations shown in Fig.@]for VT-1 (the same method 
was used for Melissa). Enough experimental data were 
taken for each source placement to minimize uncertain- 
ties from counting statistics to < 2%. 

A Monte Carlo simulation was done for each experi- 
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FIG. 5: Experimental and simulated efficiency in the 

VT-1 detector. Multiple data points at each 7-ray 

energy correspond to a different placement of the source 

in the detector. 



FIG. 4: Point source locations in VT-f (Color online). 



mentally measured source location using a point source 
and assuming isotropic photon emission. Simulated and 
experimental spectra were directly compared to deter- 
mine how well the simulation agrees with the experimen- 
tal spectrum. 

The differences between the simulated efficiencies and 
those calculated using experimental data are shown in 
Figs. [5] and [51 In VT-I, the discrepancy between the 
simulated and experimental efficiencies is significant for 
regions below the crystal. The difference increases at low 
energies. This discrepancy is likely due to uncertainties 
in detector geometries below the- crystal, which are not 
supplied by the manufacturer. The average ratio of the 
experimental to simulated efficiencies is 0.89 in this re- 
gion. The average ratio of experimental to simulated ef- 
ficiencies without considering regions below the crystal is 
0.96. For Melissa, the ratio of the experimental to sim- 
ulated efficiencies is 0.96. A constant correction factor to 
account for this difference is applied when determining 
the activity of a sample. 



C. Activity Calculations 

With the FEP detection efficiency well understood, we 
can precisely determine or set an upper limit on the ac- 
tivity of a sample. First, a sample spectrum is compared 
to the background. The background spectrum is taken 
close to the time of the assay to limit the effects of tem- 
poral variation in background, such as seasonal changes 
in radon activity. A Gaussian plus a linear background 
function is fit to each peak used for analysis. The fit is 
used to obtain the centroid of the peak and define the 
limits of the peak. The peak area is then extracted by 
summing the counts in a ±5<r width region for both the 
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FIG. 6: Experimental and simulated efficiency in the 

Melissa detector. Multiple data points at each 7-ray 

energy correspond to a different placement of the source 

in the detector. 



sample and background spectra. Subtraction of the con- 
tinuum is done using one of two methods: (I) Integrate 
the linear part of the fit extrapolated under the peak 
and subtract, (2) Define two regions, one 5er to the left 
and the other 5a to the right of the peak, average the 
counts in the two regions, and subtract it from the to- 
tal integrated peak area. If another peak is within 5a 
to the left or right, the background region with no peak 
present is used rather than the average (See Fig. [7] for 
illustration). We prefer method (2) since this does not 
rely on the goodness-of-fit. The net peak area is found 
by subtracting the background peak area from the cor- 
responding sample peak. 

The activity of a specific peak is then given by: 



A, 



net peak area [counts/second] 



(2) 



where e 7 is the peak efficiency, and m is the mass of the 
sample. 




Energy [arb] 

FIG. 7: Generic Gaussian peak as an illustration of 

integration limits. The peak is fit with a linear plus a 

Gaussian function. The limit of the Gaussian is ±5<r 

(dashed) and the linear background (solid) is 

extrapolated under the peak. 



The MaGe/Geant4 simulation determines the effi- 
ciency for detecting a decay, including geometric effects 
and branching ratios, so that the number of signal counts 
in a detector can be related to the source activity: 

signal counts = (source activity rate) x 

(counting live time) x e 7 (4) 

The sensitivity is determined by combining Eqs. ([3]) 
and ((4J, and solving for the detectable activity rate. This 
gives the sensitivity, S: 



l.Q'iy/background counts 
(counting live time) x e 7 

The sensitivity results are shown in Table IIIII 



(5) 



If no statistically significant peak is present, an upper 
limit can be placed on the activity of a sample. The 
activity at a given energy is required to be less than 
or equal to 1.6 4x ^/background counts (90% confidence 
level (C.L.)). 

Once activities for individual peaks are calculated, 
they must be related to the overall 238 U and 232 Th con- 
tent. In many cases secular equilibrium cannot be as- 
sumed, so activities or upper limits are considered sep- 
arately for 238 U (from 234 Th, 234m Pa) and 226 Ra (from 
226 Ra, 214 Pb, 214 Bi) in the 238 U decay chain and 228 Ra 
(from 228 Ac) and 228 Th (from 212 Pb, 208 T1) in the 232 Th 
decay chain. The presence of radon gas, which decays to 
214 Bi and 214 Pb, makes the measurement of 226 Ra diffi- 
cult, so measures are taken to eliminate as much radon 
as possible inside the sample cavities through the use of 
LN 2 boil-off. 



D. Estimating detection limits 

A MaGe/Geant4 simulation was performed to esti- 
mate the detection limits, defined as the net signal level 
that may be expected to lead to a detection [201, f° r a 
generic sample placed in a Marinelli beaker, modeled by a 
Teflon® ring of dimensions 14.1 cm (OD) x 8.4 cm (ID) 
x 7.1 cm (height). The detection limits for isotopes in the 
238 U and 232 Th decay chains and 40 K were calculated for 
Melissa and VT-1 using the best-achieved background 
spectrum for each detector. We define a detectable signal 
by Eq. ([3]) , where the background counts are determined 
from experimental background spectra and signal counts 
are the excess expected from the simulations of the sam- 
ple. Counts within ±5er width of the signal-peak energy 
are included. 

signal counts > 1.64 x \J background counts (3) 



V. ASSAY RESULTS 

Table IIVI shows the sample assay results to date since 
May 2008. Samples arc listed in the order that they were 
assayed. Activity limits for later assays were improved 
as a result of progress made with background reduction 
and detector performance. Disequilibrium in natural un- 
closed systems, such as plants, soil and rock, is common 
as observed in the Table Mountain rock sample. For the 
sample of aluminum stock flange coupling, disequilibrium 
was observed in the 238 U an( j 232 Th decay chains. The 
decrease in activity of the signature 226 Ra daughters have 
been observed in previous analyses of both pure and alu- 
minum alloys. This can be explained through differences 
in the chemistries of U, Th, Ra and Pb, as could occur 
through the steps required to recover aluminum from its 
ore [2l||. For lead samples, 210 Pb (T 1/2 = 22.3 • y), activ- 
ity was measured separately since it is not expected to be 
in equilibrium with the rest of the 238 U chain. The lead 
sample from Sullivan Metals had less than 2.5 Bq/kg of 
210 Pb activity. Lead bricks of unknown origin acquired 
from the University of Washington showed no measur- 
able 210 Pb signatures, but an upper limit of 10 Bq/kg 
was placed on the 210 Pb activity. 



VI. CONCLUSIONS & REMARKS 

A gamma-counting facility has been commissioned at 
KURF. The background signal rates for the Melissa and 
VT-1 detectors have been pushed to low levels. This was 
accomplished by building the facility in an underground 
location and using passive shielding, radio-pure detector 
components, and radon mitigation techniques. We have 
successfully demonstrated the analysis procedures and 
assay sensitivities required for screening materials for the 
next generation of low-background experiment. 
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TABLE III: Detector sensitivities (90% C.L.) for Melissa and VT-1 for a Teflon® ring. 







Melissa 




VT-1 


Energy 


[keV] Isotope (Chain) 


counts/day 


Detection Limit 
[mBq/kg] 


counts/dav 


Detection Limit 
[mBq/kg] 


63 


234 Th ( 238 U) 


81 ± 1 


40 


79 ± 4 


100 


93 


234 Th (238TJ) 


96 ±2 


60 


121 ± 5 


70 


1001 


234mp a (238TJ) 


2.7 ± 0.4 


50 


9 ± 1 


230 


186 


226 Ra ( 238 U) 


145 ± 3 


30 


121±5 


60 


295 


214p b (238TJ) 


117 ± 3 


8 


85 ± 4 


10 


352 


214p b (238TJ) 


114 ± 3 


5 


100 ±5 


10 


609 


214 Bi (238TJ) 


59 ± 2 


5 


53 ± 3 


8 


1120 


214 Bi (238TJ) 


13 ± 1 


8 


21±2 


30 


1764 


214 Bi (238TJ) 


9.2 ± 0.7 


8 


18±2 


30 


338 


228 Ac ( 232 Th) 


65 ± 2 


10 


48±3 


20 


911 


228 Ac (232 Th) 


5.7 ± 0.7 


3 


13±2 


10 


238 


212p b ( 232 Th) 


133 ± 3 


3 


104±5 


6 


583 


208 T1 (232 Th) 


18 ± 1 


3 


21±2 


8 


2614 


208 T1 (232 Th) 


2.9 ± 0.4 


3 


7±1 


10 


1461 


40 K ( 40 K) 


19 ± 1 


20 


32±3 


50 



TABLE IV: KURF assay results. 234 Th - activity measured in the 238 U decay series from 234 Th and 234m Pa. 226 Ra - 

activity measured in the 238 U decay series from 226 Ra, 214 Pb and 214 Bi. 228 Ra - activity measured in the 232 Th 

decay series from 228 Ac. 228 Th - activity measured in the 232 Th decay series from 212 Pb and 208 T1. Activities from 

208 T1 are divided by the branching ratio (35.94%). Not measured = NM. 



Sample 


Detector 


234 Th [Bq/kg] 


22B Ra [Bq/kg] 


228 Ra [Bq/kg] 


^ 28 Th [Bq/kg] 


4U K [Bq/kg] 


bu Co [Bq/kg] 


Table Mountain 
rock (latite) 


Melissa 


NM 


100±40 


100±40 


270±120 


790±320 


NM 


Table Mountain 
rock (latite) 


VT-1 


NM 


100±40 


100±40 


300±120 


730±290 


NM 


Supcrinsulation 
panels 


Melissa 


NM 


3.0±1.2 


NM 


0.09±0.03 


0.9±0.4 


NM 


Aluminum stock 
flange coupling 


Melissa 


7.1±2.3 


1.5±0.4 


<0.05 


1.5±0.4 


<0.2 


NM 


PMT base elec- 
tronic components 


Melissa 


<4 


1.5±1.0 


0.8±0.6 


0.6±0.4 


3.3±1.9 


NM 


PMT base elec- 
tronic components 


VT-1 


<2 


1.1±0.7 


0.8±0.4 


0.6±0.3 


3.6±1.8 


NM 


Zeolite molecular 
seive 


Melissa 


5.8±1.2 


8.2±0.8 


9.6±0.6 


10.5±0.6 


4.4±0.5 


NM 


Great Stuff™ foam 
insulation 


Melissa 


NM 


<0.2 


NM 


<0.2 


<0.3 


NM 


Axon Picocoax® 


VT-1 


<0.7 


<0.2 


0.060±0.020 


0.055±0.010 


700±200 


<0.01 


Sullivan lead bricks 


Melissa 


<0.01 


<0.002 


<0.0005 


<0.0004 


<0.003 


NM 


University of Wash- 
ington lead bricks 


Melissa 


<0.01 


<0.003 


<0.001 


<0.0004 


<0.003 


NM 


PEEK plastic 


VT-1 


<0.2 


<0.04 


<0.04 


<0.03 


<0.14 


<0.008 



